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PRECIPITATION AND WATER-YIELD RELATIONSHIPS 
ON THE LAKE GLENDALE WATERSHED 


By W. R. Boccess, R. L. RusseEtt, and A. R. GILMoRE 


Sa WATER IS A PRIMARY SOURCE OF WATER SUPPLY for many 
communities, municipalities, and cities throughout Illinois. The 
southern part of the state, with its limited ground-water resources, is 
particularly dependent upon rivers, streams, or impoundments for 
water to meet domestic and industrial needs (Illinois State Water Sur- 
vey, 19571). Except for the major rivers, stream flow is very low 
during the summer months, making it necessary to store water during 
periods of high runoff. This involves the construction and maintenance 
of dams and reservoirs. 

Siltation is a serious problem in maintaining reservoir capacity and 
may become critical in hilly sections of the state. Planning good water- 
shed management should be an integral part of reservoir planning and 
design, but this is seldom the case. Stall (1962) reported that the mean 
annual loss in reservoir storage capacity was 0.62 percent for 78 lakes 
in Illinois and represented a deposition of 1.52 tons of sediment per acre 
of watershed per year. He estimated that soil losses on these water- 
sheds could be reduced by 43 to 92 percent by the use of good manage- 
ment practices. 

Siltation can be virtually eliminated by the establishment of complete 
vegetative cover. This fact has been illustrated by two long-time studies 
in Tennessee by the Tennessee Valley Authority on the White Hollow 
(even O6L) and Pine Tree Branch (1.VA, 1962) Watersheds. White 
Hollow is in Union County and its steep topography is typical of the 
Ridge and Valley Province which lies between the Appalachian Moun- 
tains and the Cumberland Plateau. Pine Tree Branch is in Henderson 
County. Topography is hilly, and the soils are loessial silt loams inter- 
mingled with sandy coastal plain sediments of the old Mississippi Em- 
bayment. In both cases, the severly eroded, cultivated, and abandoned 
lands were reforested between 1935 and 1940. Existing woodlands were 
protected from fire and grazing, thus allowing stocking to build up to 
near normal. Sedimentation rates on both areas, with their widely 
different soils and topographic conditions, were reduced more than 90 
percent. Greatest reductions occurred during the early years following 
treatment. 


* This and similar references are listed in the bibliography. 
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The yield of water from a given drainage area is related rom 16 
amount, kind, composition, and condition of the vegetative cover. In 
most cases when forest stands are cut, increased water yields may be 
expected. Reinhart and Eschner (1962) found that in West Virginia 
hardwood stands a light cut of 1,/00 board feet per acre increased 
stream flow about 0.3 inch per year. A heavy cut of 8,500 board feet, 
however, increased stream flow ten times as much. Kovner (1956) 
found that stream flow increased 14.45 inches during the first year after 
a 39.8-acre watershed at the Coweeta Hydrologic Laboratory in the 
southern Appalachian Mountains had been clear-cut. Stream flow, 
however, declined steadily with natural regrowth of the forest cover, 
and after 13 years was only 4.99 inches greater than the stream flow 
during the base period prior to cutting. In the same area, Johnson and 
Kovner (1956) reported an annual stream-flow increase of 2.0 inches 
per year for a 6-year period following the removal of the rhododendron 
and mountain laurel understory on a 70-acre watershed. In contrast to 
these findings, there were no significant decreases in stream flow on 
the White Hollow Watershed (TVA, 1961) during the 24-year period 
following treatment despite a very great increase in the amount and 
density of the vegetative cover. IVA States?) The en cetssorean= an 
shading, reduced wind velocity, and greater humidity have evidently 
resulted in reductions in evaporation approximately equal to the in- 
creased consumptive use by plants.” (TVA, 1961, p. 44.) 

Siltation can be controlled and water purity insured by well-main- 
tained vegetative covers. The possible effects of the vegetation on 
water yields, however, should be considered in reservoir planning and 
’ design whether the over-all objective is water supply or flood control. 

Water-yield data from small watersheds, particularly those with 
well-established vegetative covers, are limited in Illinois. This is espe- 
cially true in the southern part of the state, where future growth is 
largely dependent upon development of the potentially great surface 
water resource. For this reason a study of water yields from the Lake 
Glendale Watershed was implemented in 1954 by the Department of 
Forestry of the University of Illinois in cooperation with the State 
Water Survey. The U.S. Geological Survey also participated actively 
in the project through its cooperative agreement with the State Water 
Survey in establishing and maintaining stream-gaging stations through- 
out the state. The objectives of the study were to determine annual, 
monthly, and individual storm yields as related to storm characteristics 
and soil moisture conditions. This report covers eight water years, be- 
ginning October 1, 1954, and ending September 30, 1962. 
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LAKE GLENDALE WATERSHED 


Location and Description 


The Lake Glendale Watershed (Fig. 1) is located in the SE14, S4, 
ameeuienovy 4,55, D135, R5E, Pope County, Illinois (37° 25’ N. ‘lati- 
tude; 88° 40’ W. longitude), and is adjacent to the Dixon Springs 
Agricultural Center of the University of Illinois. The area is in the 
Shawnee Hills section of the Interior Low Plateaus Province (Leighton 
et al., 1948) and is commonly known as the “Illinois Ozarks.” 

Lake Glendale has a surface area of 82 acres at spillway level and 
is maintained as a recreational area by the U.S. Forest Service. Facili- 
ties are available for camping, picnicking, swimming, and boating. 
Public use of the recreational facilities has increased greatly in the 
past decade. 

The watershed has an area of 1,350 acres (2.11 square miles), and 
about two-thirds of the land is owned by the Forest Service. Topog- 
raphy ranges from gently sloping land to steep bluffs. Seventy percent 
of the land area is on slopes less than 10 percent, 27 percent on slopes 
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Topography and location of instruments in the Lake Glendale Watershed. 
(Fig. 1) 
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between 10 and 20 percent, and 3 percent on slopes greater than 20 per- 
cent. The steepest land consists largely of bluffs along the south- 
eastern rim of the watershed. : 

In 1938, when the lake was constructed, 11 percent of the water- 
shed was cultivated, 59 percent was in unimproved pasture or aban- 
doned, and 30 percent supported stands of native hardwoods (Fig. 2). 
Erosion was moderate to severe on all lands that had been cultivated. 
The timber stands were understocked as a result of poor cutting prac- 
tices, grazing, and indiscriminate burning. 

A tree-planting program was initiated along with the lake construc- 
tion, and by 1941 forest plantations, largely shortleaf pine (Pinus 
echinata Mill.), had been established on all of the open lands owned 
by the government. The privately owned land, located largely on the 
periphery of the watershed, has either been withdrawn from cultivation 
or maintained as pasture. Land withdrawn from cultivation in this 
area soon develops a cover of broomsedge (Andropogon virginicus L.) 
which may persist for 20 years or more before it is replaced by wood 
species (Bazzaz, 1963). A broomsedge cover has high capabilities in 
preventing soil and water losses (Meginnis, 1935). The same is true 
for reasonably well-maintained pastures (Gard et al., 1959). By the 
time the present study was initiated, the entire watershed was well 
stabilized with vegetation. About three-fourths of it is now forested 
with stands of native hardwoods or plantations of shortleaf pine. 


Soils 


The soils on the watershed have developed under deciduous forest 

cover in varying depths of loess deposited over residuum or bedrock. 
Surface horizons are light-colored silt loams, and the B horizons are 
generally well developed. Permeability is moderate to very slow, and 
internal drainage ranges from poor to excessive. 

Soils that were once cultivated are primarily those of the Grants- 
burg series, with a lesser amount of Zanesville and Robbs. The Grants- 
burg soils are Gray-Brown Podzolic soils, intergrading with Red- 
Yellow Podzolic soils that developed from 80 inches of loess over 
sandstone rock. They occur on slopes ranging from 2 to 12 percent. 
The A horizon is a yellowish-brown silt loam. The B horizon is a 
silty clay loam with a subangular, blocky structure. A well-developed 
silt pan (fragipan) occurs in the lower part of the B horizon at depths 
of 30 to 36 inches (less if the soil is severely eroded). The pan is hard 
and brittle when dry, is nearly massive, and limits root penetration 
and the downward movement of water. 
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Aerial photographic mosaics of the Lake Glendale Watershed in 1938 (while 
the lake was being constructed) and in 1959, (Fig. 2) 
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The associated Robbs series are imperfectly drained Planasols inter- 
grading with Gray-Brown Podzolic soils that occur on slopes of 1 to 4 
percent. The silt loam A horizon may be up to 16 inches thick. The 
silty clay loam B horizon has a higher clay content than that of the 
Grantsburg and is slowly permeable to water. Because of the low slope 
gradient, erosion has never been a serious problem with these soils. 

The Zanesville series developed in about 40 inches of loess over 
sandstone. Profile development is somewhat similar to that of the 
Grantsburg soils, but is shallower, and the fragipan development is 
weaker or entirely missing. 

Practically all of the better soils were cleared for agriculture during 
the early days of settlement. Only those that were either too steep, 
rocky, or shallow were excluded from cultivation. Thus soils under the 
native hardwood stands are in this category. They largely belong to 
the Wellston and Muskingum series, a complex of Muskingum and 
Wellston, and the steeper phases of Zanesville. There is also some 
rough, stony land and rock outcrops, but these are more of a “land 
class” than a definite soil series. 

Soil characteristics that are hydrologically important include (1) the 
slowly permeable fragipan of the Grantsburg soils, (2) the shallow 
profiles of the soils occupied by native timber, and (3) the amount of 
erosion of the formerly cultivated soils. In general, these factors 
have resulted in soil profiles that have relatively low moisture-storage 
capacities. Thus they contribute heavily to the precipitation-runoff 
characteristics of the watershed. 


COLLECTION OF FIELD DATA 


A network of five standard and three recording rain gages, well 
distributed over the watershed, was used to measure precipitation 
(Fig. 1). However, only four standard gages were in operation for the - 
entire 8-year period; data from these were used in the analyses. The 
recording gages were equipped with 12-hour charts and provided in- 
formation on storm intensity and duration. Data were collected as soon 
after each storm as practicable. 

Precipitation was measured in the open at each gage site. No allow- 
ance was made for interception by the vegetation, even though such 
data are available for shortleaf pine plantations in this area (Boggess, 
1956). This procedure was followed because all data from weather 
bureau stations are based on open gages, thus making possible the use 
of the precipitation-water yield relationships of this watershed in other 
similar locations. 
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Two gaging stations were established to determine flow into the 
lake and to record water levels. These installations were made by the 
U.S. Geological Survey and the State Water Survey as a part of their 
statewide stream-gaging network. Each station consisted of a Stevens 
Stage Level Recorder, installed over a stilling well in a USGS standard, 
corrugated iron shelter. 

The headwater station was specifically designed to measure low 
flows over a concrete control structure. The drainage area for this 
station was 1.04 square miles, or approximately half of the entire 
watershed. 

The lake-level recording station was located near the spillway and 
provided continuous records of water levels to the nearest 0.01 foot. 
The instrument was not effective when the water levels fell to more 
than 3 feet below the spillway. 

Heights recorded by the instruments were checked each week by 
visually reading staff gages that had been established at each station 
and referenced to spillway level. At the same time each instrument 
was checked for possible malfunction and the observation time was 
noted on the strip chart of the stage recorders. This permitted adjust- 
ment of time over the month-long period represented by each chart. 

Rating tables to translate gage heights into water discharge were 
prepared by the USGS from stream-gaging data collected over a wide 
range of water levels at each station. However, these tables could not 
be used to determine discharge from the entire watershed when the lake 
was below spillway level. For these conditions discharge was calculated 
from equivalent gage height-lake volume curves prepared from surveys 
of the lake basin topography. The underwater elevations were deter- 
mined from soundings along permanently marked range lines. The 
420-, 422-, and 424-foot contour lines were established by the Topo- 
graphic Division, U.S. Geological Survey, and allowed the extension 
of the capacity tables to water levels above that of the spillway. The 
spillway level was 418.38 feet above sea level. 

Data from the lake-level gage are incomplete for the water years 
of 1957-58, 1959-60, 1960-61, and 1961-62, when the lake was lowered 
during the winter months by the U.S. Forest Service to repair facilities. 
Storm yields during these periods were estimated by predicting equa- 
tions which will be discussed later. 

Ground-water levels were measured in two wells, one located in a 
pine plantation and the other in an abandoned field (Fig. 1). The wells 
were dug to bedrock (10 to 12 feet) and cased with 10-inch bell tile. 
Continuous-stage level recorders were placed on each well. The re- 
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corded heights were checked by visual measurements of water levels in 
the well. 

Gravimetric soil-moisture determinations were made weekly adja- 
cent to each ground-water well. These were discontinued after two 
years, as it was possible to use the ground-water levels as an index of 
soil moisture storage capacity. 


ANALYSIS OF DATA 


Precipitation 


The precipitation data collected from the four standard gages were 
tested by a double-mass curve technique to find whether any inconsis- 
tencies existed in the record of any one gage (Schneider and Ayer, 
1961). The results are shown in Fig. 3. This curve was constructed by 
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CUMULATIVE PRECIPITATION (INCHES) — AVERAGE 4 STATIONS 


Relationship between cumulative precipitation for gages in the study area 
and cumulative precipitation from average of four stations. (Fig. 3) 
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plotting the cumulative precipitation for one gage against the average 
for the group of gages during the same period of record. Since the 
plotted curves showed no abrupt change in slope, it was assumed that 
there were no inconsistencies in the data. The unweighted mean of the 
four gages was used to express the equivalent uniform depth of pre- 
cipitation over the watershed. The consistently low readings of the 
gage on the spillway were probably related to its greater exposure to 
wind. 

In the 8-year period of record for the watershed, 727 storms (mea- 
surable precipitation occurring 6 hours or longer after the preceding 
storm) were recorded. This is an average of 91 storms per year. The 
number of storms occurring in a given water year varied from 66 
storms in 1955-56 to 108 storms in 1954-55 (Table 1). 

The average annual precipitation on the watershed was 48.66 inches 
for the 8-year period. This is slightly above the 24-year recorded aver- 
age of 46.64 inches at the Dixon Springs Agricultural Center, located 
about 1 mile north of the watershed. Annual precipitation varied from 
37.75 inches to 68.30 inches during the period of record on the water- 
shed. These amounts were within a 5-year recurrence interval for all 
years except 1957-58, when the total was 68.30 inches. This amount 
would be expected to occur only once in 25 years (State Water Survey, 
1957). The precipitation data are summarized in Table 2. 

The minimum monthly total of 0.46 inch occurred in December, 
1955, and the maximum of 16.30 inches in July, 1958. Monthly precip- 
itation totals of less than 1 inch are not uncommon in this area, but the 
maximum total of 16.30 inches has been exceeded only twice in the 
history of weather bureau records for this area in southern Illinois. 


Table 1.— Yield-Producing Storms and Evapotranspiration and Losses 


Number Percentage Evapotrans- 
Number of storms of storms piration 

pieper year of storms producing producing plus losses 
yield yield (inches) 
ao >. 5 rr 108 53 49 3032 
LS SS a 66 35 53 27.39 
DS 100 56 56 33305 
US) Se 105 63 60 41.09 
ee IN oe nes ae cine eo 2 80 44 55 33470 
LOCC) 71 39 55 26.40 
CC 43 ne 100 50 50 32.86 
SOG, 97 36 il 30.50 
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New Burnside and Golconda reported 16.94 and 18.52 inches, respec- 
tively, in January, 1937, and 17.59 inches was recorded at the Dixon 
Springs Agricultural Center in January, 1950. The recurrence interval 
for the July, 1958, total is probably once in 75 to 100 years. 

The average storm size was 0.54 inch, ranging from .01 (trace) 
inch to a maximum of 5.08 inches (July, 1958). This maximum storm 
has a recurrence interval of once in 18 years (Huff and Neiteio59 
Storm size and the distribution of storms by day of occurrence are 
shown in Figs. 9-16, pages 27-34. The distribution and size of the indi- 
vidual storms varied considerably, although yearly and monthly totals 
compared favorably with long-term averages. Storm variation during 
the 8-year period, however, is perhaps not abnormal. 

Intensities for the 727 storms ranged up to 3.30 inches per hour, 
which occurred during the storm of July 20, 1955. Several storms had 
higher intensity rates but for shorter intervals of time. 


Reduction of Hydrographs 


Since an objective of this study was to determine individual storm 
yield, the normal method of hydrograph reduction, in which daily 
yields are calculated (U.S. Department of Agriculture, 1962), could 
not be used because of carryover of storm yield from one day to the 
next. Therefore each storm was analyzed separately to determine 
individual storm yield. 

Yield of each storm was reflected by a rise in the lake level. The 
rise for each storm was corrected for precipitation on the surface of 
the lake by subtracting the volume of water as read from the gage 
height-capacity curve that was equivalent to the amount of precipita- 
tion. Thus the yield from the watershed, or land area, could be deter- 
mined. 

Storm yields fell into four classes: (1) storms in which the lake 
rise did not exceed spillway level; (2) storms during which spillway 
level was exceeded; (3) storms when lake levels were already above 
spillway level; and (4) storms whose discharge overlapped that of pre- 
ceding storms. 

The yield of storms in which the lake rise did not exceed spillway 
level was computed from gage height-volume curves to the nearest 
thousand cubic feet of water. 

Storm yields during which spillway level was exceeded were com- 
puted by a combination of the gage height-capacity relationship and 
the integration of the storm hydrograph. This procedure is outlined 
below. 
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® Volume of lake rise from start of storm to spillway level was 
computed from the gage height-volume curves. 

@ Discharge over the spillway was calculated by first integrating 
the storm hydrograph into segments of one-half hour up to 20 hours 
duration, the length of time depending on changes in the hydrograph 
slope. The average gage height during each segment was determined, 
and discharge in cubic feet per second for this height was obtained 
from the rating table based on gage height-discharge relationship. The 
cubic-foot-per-second discharge was then corrected to a cubic foot 
volume figure. This second step was required only when the lake was 
above spillway level at the beginning of a storm. 

® Total storm yield was equal to the volume increase up to spill- 
way level, plus the discharge over the spillway. 


In calculating the yields from overlapping storms, it was first meces= 
sary to separate flow from the first storm from that of the second. This 
involved approximating the recession leg of the first storm on the basis 
of that actually recorded for the combined storms. This was done by 
beginning at the point on the recession leg of the first storm where 
runoff from the second storm first appeared, and then constructing the 
remainder of the first recession leg parallel to that shown on the hydro- 
graph for the combined storms. The method was simple but had a 
tendency to overestimate the yield from the first storm, particularly 
during periods when the soil was wet and subsurface flow was high. 
This was probably related to the fact that the recession leg on the hydro- 
graph showed the effect of both storms, particularly on subsurface 
flow (along the top of the silt pan). Thus the estimated recession leg 
of the first storm was lengthened somewhat, because it was drawn 
parallel to that of both storms combined. The overestimation of yield, 
however, was not great for an individual storm and was compensated 
for in total yields by a like reduction in the yield for the last storm in 
an overlapping series. A procedure somewhat similar to the one de-- 
scribed above was reported by Reinhart (1964) but appeared after our 
analyses were completed. 


Regression Analysis 

Equations were developed to express water yields in relation to 
precipitation. However, some of the storms that fell during the 8-year 
period could not be utilized in the analyses. These included 108 storms 
in the 1954-55 water year, when ground-water records were not avail- 
able, and 99 storms that occurred when the lake was drained. After 
considering several factors including storm intensities, antecedent pre- 
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cipitation, days since last significant storm, and ground-water levels, 
it was found that ground-water level was the most consistent and reli- 
able single factor that could be used to predict whether or not water 
yield might be expected from a given storm. 

Data from the ground-water well in the southwest corner of the 
watershed (Fig. 1) were used because its location more nearly repre- 
sented the watershed, both from the standpoint of soils and of vegeta- 
tion, than did the well in the abandoned field near the headwater. The 
specific ground-water level used in connection with a given storm was 
the low point on the recorder chart prior to the rise resulting from sub- 
surface storm flow. 

When ground-water levels were plotted against storm yields, it 
appeared that yield was related to levels either above or below a point 
approximately 4 feet below the ground surface (Fig. 4). Water yields 
occurred for most storms greater than 0.20 inch in size when ground- 
water levels were less than 4 feet from the surface. At levels greater 
than 4 feet, yield occurred only when the storage capacity of the soil 
was exceeded during the storm, or when storm intensity was greater 
than the infiltration rate. 

The 4-foot ground-water level was an index of soil moisture storage 
capacity. For the Grantsburg soil, the most abundant on the water- 
shed, it indicated the point where the soil horizons were wet to field 
capacity or above, down to the fragipan. Additional precipitation either 
(1) ran off the surface (overland flow), (2) percolated downward to 
the fragipan and moved down the slope along the top of this slowly 
permeable layer, or (3) slowly wet the pan itself and the rest of the 
soil mantle down to bedrock. A rise in ground-water levels was largely 
dueto1(Z): 

From the standpoint of water-yield prediction, the storms were 
grouped into three classes: 


1. Storms falling on “wet soil,” with no moisture storage capacity 
available in the profiles (ground-water levels less than 4 feet from the 
surface). There were 330 storms in this category, and 191 resulted in 
water yields (Area 1, Fig. 4). 

2. Storms falling on “dry soil,’ with moisture storage capacity avail- 
able as indicated by ground-water levels greater than 4 feet from the 
surface. This group included 279 storms; those having some yield are 
shown in Area 2, Fig. 4. 

3. Transitional storms that began on “dry soil’ but during which 
the moisture storage capacity of the soil was exceeded. Only 10 storms 
fell in this category (Area 3, Fig. 4). 
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DEPTH TO GROUND WATER (FEET) 


Relationship between water yield and depth to ground water. 
(Fig. 4) 
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Water yields for the various storms were plotted against storm 
size and a linear regression line was fitted to these points, according to 
the following formula: 


Y =y-+4+ b (X —&), where 


Y = point on the regression line or the predicted value of the water 
yield, 

y = mean of water yield, 

x = mean of precipitation, 

X = amount of precipitation (uniform depth over entire water- 


shed), 
b = the regression coefficient. 


In each case an analysis of variance was used to test the significance 
of the regression relationship. 

The regression equation Y = 0.91X — 0.16 defined the precipita- 
tion-water yield relationship on “wet soil” (Fig. 5). The equation was 


WATER YIELD (INCHES) 


00 8 2 ‘ ¥ = O0:91Xu= OU6 
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Relationship between water yield and precipitation for storms on wet soil. 
(Fig. 5) 
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significant at the 1-percent level. The correlation coefficient (r) was 
0.95 and the coefficient of determination (1?) was 0.90. 

The predicting regression equations for a precipitation-water yield 
relationship often yield a result of zero even though precipitation is not 
zero. This occurs because it takes a certain amount of precipitation to 
exceed the moisture storage capacity of the vegetation and allow 
throughfall to reach the ground. 

Water yields resulted from about one-fourth of the storms that fell 
on “dry soil,” and these yields were relatively small. The relationship 
between water yield and precipitation for the “dry soil” condition was 


Y = 0.035X — 0.03. 
While the relationship was significant, the correlation coefficient (0.50) 
was much weaker than that for the “wet soil” storms. Thus the in- 
crease in yield was small in comparison with increase in precipitation. 
The transition storms were first classified as “dry soil” storms, since 
soil moisture storage was available at the beginning of the storm. How- 
ever, this caused an overestimation of yield from storms falling on 


Y =051X-8.40 


WATER YEAR RUNOFF (INCHES) 


e@= Partially Estimated 


O 
25 30 35 40 45. 50 55 CURR lma” 
WATER YEAR PRECIPITATION (INCHES) 


Relationship between water yield and precipitation by water year for eight 
water years beginning October 1, 1954. (Fig. 6) 
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dry soil. As there were only 10 transition storms, compared with 279 
“dry soil” storms, it seemed advisable to develop an equation for the 
transitions. The relationship could be expressed as 


A 


Y = 0.416X — 0.069. 


It was significant at the 5-percent level and had a correlation coefficient 
On10 72, 

The equation expressing the precipitation-water yield relationship 
for the eight water years (Fig. 6) was 


Y = 0.51X — 8.40, 


which was significant at the l-percent level. The correlation coefficient 
(r) was 0.93 and the coefficient of determination (r?) was 0.86. This 
equation could be used to predict water field from annual precipitation 
ranging from 37.75 to 68.30. These amounts cover a recurrence interval 
of 5 years for the minimum and 25 years for the maximum yearly 
amounts. 


WATER YIELD 


Of the 727 storms measured during the 8-year period, 376, or 52 
percent of the total number, produced some yield from the watershed. 
The number of storms producing yield varied from 35 to 63, with an 
average of 47 per year (Table 1). 

Water yield was composed of both surface (overland) and sub- 
surface flow. No attempt was made to separate these two elements, 
since total yield was the main objective of the study. However, sub- 
surface flow was a major contributor to total yield, due both to the 
presence of the slowly permeable fragipan and the relatively low mois- 
ture storage capacities of the soil profiles. 

Subsurface flow was believed to be largely made up of down-slope 
seepage along the top of the fragipan, since a perched water table 
formed there during prolonged periods of wet weather. Deep seepage 
(moisture that moves downward to bedrock and then into drainage 
channels) may have made some contributions to subsurface flow but 
the amount was probably small in comparison with the movement of 
perched water along the pan. 

A summary of the monthy and annual water yields from the water- 
shed is shown in Table 3. The totals are based on individual storm 
yields which eliminated monthly carryover, since the yield of each 
storm was recorded in the month during which it began. 

Annual water yields varied from 9.10 to 27.21 inches, with the 
lowest yield in 1954-55 and the highest yield in 1957-58. The average 
was 16.74 inches, or 34 percent of the average annual precipitation. 
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Monthly water yields varied from none to 8.27 inches. There were 
no water yields largely during the fall and early winter months when 
the soil profiles were in the process of recharge following depletion 
during the growing season by evapotranspiration. The maximum 
amount occurred in July, 1958, when precipitation was far in excess of 
normal. At the beginning of the month, ground-water levels were 
following a normal downward trend, which was reversed with storms 
that produced 3.3 inches of rainfall on July 6 and 7. From July 15 to 
July 21 more than 12 inches of rain occurred in storms from 0.10 to 
5.08 inches in size, and ground-water levels reached the highest point 
recorded during the 8-year period. Soil profiles were completely water- 
logged, and about 50 percent of the precipitation was returned to the 
lake as surface or subsurface runoff. Average July yields for the rest of 
the study period amounted to 7 percent of the precipitation. 

The seasonal aspects of water yield are shown in Fig. 17, page 35. 
Monthly averages expressed as a percentage of precipitation are shown 
in Table 1. Winter and spring had high yields and may be considered a 
period of moisture surplus. Normally the low-yield period occurred in 
the summer and fall and was composed of two distinct periods. The 
first was a period of depletion that resulted when evapotranspiration re- 
quirements exceeded precipitation and moisture was withdrawn from 
soil storage. It was closely related to the growing season. The second 
period was that of recharge. It occurred at the end of the growing 
season and continued until the soil moisture storage capacity was Satis- 
fied. The beginning of the recharge cycle usually marks the low point 
of soil moisture storage and stream flow during the year. For this 
reason October 1 was selected as the beginning of the water year. This 
date is generally used as a basis for compiling hydrologic data in the 
United States. 

The annual soil moisture trends are illustrated by the plot of ground- 
water levels in Figs. 9-16. The periods of moisture excess, depletion, 
and recharge are easily recognized. The moisture excess period was 
marked by water levels less than 4 feet from the ground surface, deple- 
tion by a downward trend, and recharge by an upward trend of the 
curve. The beginning and ending of each period varied considerably 
from year to year, depending largely on precipitation patterns. All 
three periods — excess, depletion, and recharge — occurred in July, 
1958 (Fig. 12, page 30). July is normally characterized by rapid mois- 
ture depletion. 

The relationships between storm size, ground-water levels, and water 
yields are also shown in Figs. 9-16. Yield from a given storm depended 
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largely on whether precipitation fell on dry soil or wet soil as previously 
discussed in relation to height of the ground-water table. Occasional 
storm flow occurred when the soil was dry. This resulted from storms 
of very high intensity that caused surface runoff, or when the soil-mois- 
ture storage capacity was reached during a storm. 

The relationships of monthly precipitation, water yield, and ground- 
water levels for the eight years of record are shown in Fig. 7. 


Relationships Between Hydrologic Factors 


Mass curves were used to study the relationships between hydrologic 
factors on the watershed (Fig. 7). When accumulated values for an- 
nual precipitation, water yield, and evapotranspiration plus losses are 


ACCUMULATED HYDROLOGIC FACTORS (INCHES) 
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Relationship between accumulated hydrologic factors and water year. 
(Fig. 7) 
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plotted against time, the slope of the curves indicates the rate of accu- 
mulation of each factor, and changes can be more readily detected than 
from tabular values (TVA, 1963). 

The high precipitation in 1957-58 was reflected by an increase in 
slope, while a corresponding decrease resulted from the below normal 
precipitation in 1959-60. 

The yield curve has three different slopes, each lasting about three 
years. The greatest increase in 1957-58 resulted from the high yields 
for July, 1958. Other changes were mainly due to seasonal distribution 
of precipitation. 

Evapotranspiration and other losses represent the difference between 
precipitation and water yield. Included are interception by the vegeta- 
tion, evaporation from the ground surface, transpiration, and losses due 
to deep seepage (moisture used in wetting the soil from the fragipan to 
bedrock). The greatest slope change is again related to the wet year of 
1957-58, when more water was available for evapotranspiration than in 
other years. The slope for the last three years is about the same as that 
at the beginning of the study. Maximum losses were 41.09 inches in 
1957-58, and the minimum of 26.40 inches occurred in 1959-60 (Table 
1). Average annual losses from the watershed were 31.92 inches. 
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factors. , (Fig. 8) 
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Double-mass curves were also plotted for water yields and evapo- 
transpiration plus losses to determine any possible changes in the rela- 
tionships between the hydrologic factors for the period of study (Fig. 
8). The slope of the evapotranspiration-plus losses curve is about con- 
stant, although a small change occurred in 1957-58. More changes, 
however, are evident on the water-yield curve. This is related to the 
fact that the ratio between annual precipitation and water yield is not 
strictly linear, as the ratios appear to be higher in wet years than in dry 
years. The slope of the double-mass curve might also be affected by 
changes in the land-use pattern or in the vegetation. However, these 
two factors appeared to be relatively unimportant. Land-use patterns 
are known, and the vegetation cover on the watershed was well estab- 
lished at the beginning of the study. While density of the vegetation no 
doubt increased due to normal growth, it did not cause an apparent 
increase in evapotranspiration. At least such changes are not apparent 
on the curve; thus variations in the slope of the water-yield curve are 
probably related to the amount of precipitation and its seasonal distri- 
bution during each year. 


DISCUSSION 


The results of this study emphasize some important points in the 
management of watersheds. First of all, yields from this forested 
watershed are essentially the same as those given for the general area 
in the Atlas of Illinois Agriculture (1958), which were based on 20- to 
40-year records from a number of stations. The average annual yield 
_ of 16.74 inches is equivalent to approximately 455,000 gallons per acre. 
With an average daily use of 50 gallons per capita, each acre would 
supply about 25 persons, or the 1,350 acres would support a population 
of about 34,000 people. Of paramount importance is the fact that the 
water in Lake Glendale is usually clear and seldom more than murky. 
Little treatment, other than chlorination, would be required to make it 
potable. 

Based on observations made when the lake was drained, siltation 
has been negligible during the 25-year life of the reservoir. Low-cut 
stumps in the basin were clearly visible and showed little evidence of silt 
accumulation. Thus the reservoir should have a long period of 
usefulness. 

Adequate water yield, clear water, and control of siltation are the 
main objectives of watershed management. The latter two are a direct 
result of the excellent cover of vegetation. The yields were obtained in 
spite of the vegetation which transpired large quantities of water. 
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As the potential surface water resources of southern Illinois are 
impounded for domestic use, it is hoped that the reservoir planners and 
developers will be cognizant of the values that accrue from establishing 
and maintaining good covers of vegetation on the watersheds. The kind 
of vegetation, either grassland or forest, is relatively unimportant in this 
area of relatively shallow soils. Either will prevent erosion, provide 
clear water, and improve the general infiltration capacities of the soil. 
The salient point is to make reservoir development and watershed man- 
agement a joint operation. 


SUMMARY 


Hydrologic factors on the Lake Glendale Watershed were measured 
for eight water years, beginning October 1, 1954, and ending September 
30, 1962. This 1,350-acre watershed was well covered by vegetation, 
with forest making up about two-thirds of the total area at the beginning 
of the study. Little vegetation change, other than normal growth, oc- 
curred during the measurement period. 

The hydrologic data were analyzed to determine annual, monthly, 
and individual storm yields as related to storm characteristics and soil 
moisture conditions. 

In the period of record, 727 storms were recorded on the watershed 
with an average of 91 storms per year. The average annual precipita- 
tion was 48.66 inches, which was 2.02 inches above the long-time aver- 
age. Monthly precipitation varied in distribution from year to year with 
the range of 0.46 to 16.30 inches. Average monthly totals agreed closely 
with the long-time averages. The precipitation of 16.30 inches in July, 
1958, was very unusual and has a recurrence interval of approximately 
75 to 100 years. The maximum storm size was 5.08 inches, which 
occurred during July, 1958. Hourly storm intensities varied from a 
trace to 3.30 inches. 

Each storm was analyzed separately to determine individual storm 
yield. Yields were reflected by a rise in the lake level due to flow from 
the watershed and precipitation on the lake surface. A method was 
devised to estimate storm flow when storm yields overlapped in time. 
The equivalent volume of water was computed by stage height-volume 
tables and stage height-discharge tables. The volume was corrected for 
the precipitation on the lake surface to find water yield from the land 
area only. 

Water yield to the lake was made up of three component flows: 


1. Overland flow, when infiltration capacity was exceeded by the 
intensity of precipitation and surface flow occurs. 
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2. Perched flow, when the percolating water reached the fragipan, 
much of it flowing along the pan to be contributed to the lake. 

3. Base flow (from water below the fragipan). When the pan was 
wetted, downward percolation occurred below the fragipan and the 
source of base flow from ground water was recharged. 


Ground-water levels were found to be an index of soil moisture or 
storage capacity. When the ground-water levels were less than 4 feet 
from the surface (‘wet soil’’), the soil had little or no moisture storage 
capacity. When the ground-water levels were more than 4 feet from 
the surface (“dry soil”), storage capacity was available and moisture 
was required to recharge the profile. Ground-water levels that defined 
soil moisture conditions were the most reliable single factor studied that 
could be used to predict water yield. 

Storms for water-yield prediction were separated into three pre- 
cipitation-water yield relationships: 

1. Storms falling on “wet soil,” represented by the regression equa- 
tion Y = 0.91X — 0.16. 

2. Storms falling on “dry soil,” represented by the equation Noi 
O35 Stee U3! 

3. Transition storms during which the soil moisture storage capacity 
was exceeded, represented by the equation Y = 416X — .069. 

A uniform depth of precipitation over the entire watershed was as- 
sumed in each case. 


A water-year precipitation-water yield relationship was developed. 
This relationship could be represented by the equation Y =.51X — 8.40, 
- which was tested statistically and found significant at the 1-percent 
level. 

Measureable yield was produced from 376 storms or 52 percent of 
the total number of storms. Water-year yield-producing storms varied 
from 35 to 63, with an average of 47 per year. The water-year totals. 
varied from 9.10 to 27.21 inches, with an average of 16.74 inches. 
Monthly water yields varied from none to 8.27 inches. Water yields 
were seasonal with the high yields occurring in the winter and spring 
and the low yields occurring in the summer and fall. 

Changes in the precipitation-water yield relationships for the period 
of study were related to changes in amount and distribution of precipi- 
tation during each year. 
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